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ABSTRACT 

We study the feedback from an AGN on stellar formation within its host galaxy, 
mainly using one high resolution numerical simulation of the jet propagation within 
the interstellar medium of an early-type galaxy. In particular, we show that in a re- 
alistic simulation where the jet propagates into a two-phase ISM, star formation can 
initially be slightly enhanced and then, on timescales of few million years, rapidly 
quenched, as a consequence both of the high temperatures attained and of the reduc- 
tion of cloud mass (mainly due to Kelvin- Hclmholtz instabilities). We then introduce 
a model of (prevalently) negative AGN feedback, where an exponentially declining 
star formation is qu enched, on a very short tim e scale, at a time tAGN, due to AGN 
feedback. Using the iBruzual fc Chariot! (|2003f ) population synthesis model and our 
star formation history, we predict galaxy colours from this model and match them to 
a sample of nearby early-type galaxies showing signs of recent episodes of star forma- 
tion (Kaviraj et al. 2007). We find that the quantity t ga i — tAGN-, where t ga i is the 
galaxy age, is an excellent indicator of the presence of feedback processes, and peaks 
significantly around t ga i — tAGN ~ 0.85 Gyr for our sample, consistent with feedback 
from recent energy injection by AGNs in relatively bright (Mb <i —19) and massive 
nearby early- type galaxies. Galaxies that have experienced this recent feedback show 
an enhancement of 3 magnitudes in NUV(GALEX) — g, with respect to the unper- 
turbed, no-feedback evolution. Hence they can be easily identified in large combined 
near UV-optical surveys. 

Key words: galaxies: jets - intergalactic medium - galaxies : elliptical and lenticular, 
cD - galaxies : evolution. 



1 INTRODUCTION 

Active galactic nuclei (AGN) have been advocated in re- 
cent years as sources able to influence the evolutionary 
history of stellar populations within their host galaxies. 
They are ubiquitous and lie in the cores of a wide vari- 
ety of galaxies, possibly suppressing, or possibly on occa- 
sion speeding up, star formation (SF). Massive and bright 
early-type galaxies (ETGs) have been shown to host very 
old stellar populations where SF stopped early, while fainter 
ETGs show lower ages and more protra cted SF phases 
(|Thomas et al.ll2005l . lDe Lucia et aLlhoOfj ). This is the so- 
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called downsizing scenario, which is predicted by observa- 
tions and generated in s o me recent simulations (ICowie et al 
I l996| . iBorch et alJlioOol iBundv et al.ll2006l. |Pe Lucia. et"al. 
2009. Uragerl l200ol iThomas et al.ll2005l . iNelan et all 120051 
iTortora et al.l 120091 ). However, only very recently semi- 
analytic simulations have been able to reproduce some of 
the k ey observations l|Cattaneo et al.l [2006. Dc Lucia et alj 
l200rJ V 



iKauffmann fc Chariot! l|l998h predicted results that 
contrast with the downsizing scenario. The main reason for 
this discrepancy may be that supernova feed back in mas- 
sive g alaxies is not suffi cient to quench SF (|Pekel &: Silld 
1 19861 . iBenson et all 120031 ). Feedback from AGN is the ad- 
ditional source able to heat the cold gas that fuels SF in 
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galaxies. It is such effects that allow us to match the evolu- 
tion of massive galaxies predi cted from simulations with re- 
sults from observations (e.g.. [Kavirai et al.ll2005l. ISilkll2005l . 



iBower et aDl200d . lCroton et al.ll2006l , lDe Lucia et al.ll200d) 



However, these studies did not account for the downsizing 
phenomenon. 

Feedback from AGN plays a major role during differ- 
ent phases of galaxy evolution. Higher activity is observed 
at high redshifts. At these early phases of galaxy evolution, 
during a major merger, the gas accreted by central super- 
massive bla ck holes is ejected within the host galaxy, a nd 
inhibits SF i|Di Matteo et al.ll2005l . ISpringel et alj|2005al lbh. 
At more recent epochs, the effect of AGN feedback is 
weaker, but still su fficient to quench SF, particularly near 
the central regions (|Croton et all 120061 . IBower et al.l I2006L 
ISchawinski et al.ll2006l . 2008). Evidence for the presence of 
AGNs is evident i n ma ny previous studies. In particular, 
Kauffmann et al. 1 (|2003h collected more than 20000 narrow- 
line AGNs from SDSS at z £ 0.3. Type-II AGN are found 
within host galaxies having structural properties similar 
to normal ETGs, and reside almost exclusively in massive 
galaxies with M+ > 10 10 Mq . Many of these galaxies have 
experienced a recent burst of SF or show evidence of ongoing 
SF, thus indicating that a black hole and prominent source 
of fuel supply are necessary ingredients for AGN feedback, 
but n o evidence of SF quenching is found (|Kauffmann et al~l 
While the existence of a black hole within massively 
star forming galaxies is a rarity in our cosmological neigh- 
bors, the situation is different at higher redshifts, w here 
the AGN frequency is higher l|Hasinger et alj 120051 ) and 
the feedback could be effective in producing red galaxies 
jBenson et al.ll2003l , ICattaneo etaflbood , ISchawinski et al.l 
l200Sh . 

More recently, a series of analyses have been devoted 
to the detection of recent star formation (RSF) within the 
last Gyrs of galaxy lifetimes. Ultraviolet data allowed the 
detection of small fractions of young stars, allowing one 
to tra ce this residual SF in ETGs both in the local uni- 
verse llYi et al.1 2005, Ka yirai et al. I l2007al ) and at high red- 
shift (|Kavirai et al.l(200Sf ). The link between SF and AGNs 
is strengthened by other analyses, where AGNs have been 
found to be responsible for the migration of star-forming 
galaxies, lying within the bl ue cloud, to almost qu iescent 
locat ed in the red seque nce (ISchawinski et all [20 06. 2007, 
2008. iBildfell et alj|200g| ). In lKavirai et al.1 (l2007bi Tthe re- 
cent feedback effect in E+A galaxies is analyzed and the 
bimodal correlation of the quenching efficiency with mass 
and luminosity is interpreted as due to negative feedback 
from supernovae or AGNs. 

The physics of the SF-AGN connection is still poorly 
understood, and the empirical models of SF quench i ng by 
AGN used in diffe rent works l|Granato et al.1 l200ll . |2004 
iMartin et alj [20 07a). although reasonable, do not yet have 
a solid motivation. Thus, a primary goal should be to 
analyze this phenomenology by means of hydrodynamical 
simulations, modelling the propagati on of jets produced 
by A GNs within a gaseous medium i|Scheuerl fl97i . iFalld 
Il99ll ) and their i nteraction with an inhomogeneous in- 
terste l lar medium dSaxton et al. 20051 . iKrause fc Alexander! 
12001 Sutherland & Bick nelll 120071 )" In a previous paper 
( Antonuccio-Delogu fc Silkll2008l . paper I hereafter) we used 
an Adaptive Mesh Refinement (AMR) code to follow the 



evolution of the cocoon produced by the jet propagating in 
the ISM/IGM: we have analyzed the thermodynamic evolu- 
tion of a cloud embedded within the cocoon and seen how 
its SF is modified. We found the SF to be quenched on time- 
scales of 10 5 - 10 6 yrs. 

In the present paper, we deduce a Star Formation His- 
tory (hereafter SFH) directly from a numerical simulation 
of the interaction of an AGN with the multiphase ISM of 
its host galaxy. In principle, a realistic simulation should at 
least encompass a range of scales spanning more than 5 or- 
ders of magnitudes in the spatial coordinates, down to single 
star-formation scales, and should include a very wide range 
of heating and cooling phenomena. In this work we have re- 
stricted ourselves to consider only the mechanical feedback 
from the AGN, and its effect on the thermodynamic state 
of the ISM, thus neglecting the radiative feedback. We have 
modelled the ISM as a two-phase, multicloud system, and we 
have followed the thermodynamic evolution of both phases, 
particularly of the cold, star-forming phase. We introduce 
a feedback model where negative feedback, i.e. the suppres- 
sion of SF within the volume affected by the jet, plays a 
dominant role. Using an empirical prescription to take into 
account the effect of cocoon on SF within cold clouds, we 
present a scenario where SF in unperturbed ETGs is de- 
scribed by an exponential SF law, corresponding to an early 
starburst followed by a slow decline. This evolution, starting 
at some time tAGN is quenched by a feedback process, and 
we follow the evolution of galaxy colors, determining the free 
parameters such as galaxy age t aa i and tAGN by comp aring 
with a previously studied sample l|Kavirai e t al. 2007j]|). 

For the first time, the amplitude and typical timescales 
of negative AGN feedback are directly derived from a simu- 
lation designed to address this phenomenology. 

The plan of the paper is as follows. In Sect. [2] we dis- 
cuss the simulation set-up, while in Sect. Owe describe the 
propagation of the cocoon and its influence on SF. In order 
to compare predicted results wit h observations, we use sin - 
gle burst population models from lBruzual fc Chariot! (2003), 
suitably convolved with the SF rate predicted by our sim- 
ulation: we describe our procedure and predicted synthetic 
colours in Sect. 3] Finally, Sect.[S]is devoted to a comparison 
with observations, and conclusions are detailed in Sect. [6] 

Whenever required, we will adopt the 3-year WMAP 
"concordance" model: h=0.74, Q, m = 0.234, Q b h 2 = 0.0223, 



corresponding to a Universe age of t\j n 



13.6 Gyr. 



2 SIMULATION SETUP 

To perform the simulation, we used FLASH v. 2. 5 
(|Frvxell et alj |2000). a parallel, Adaptive Mesh Refinement 
(AMR) code, which implements a second order, shock- 
capturing PPM solver. The modular structure of FLASH 
allows the inclusion of various physical effects, includ- 
ing a source of external heating, radiative cooling and 
thermal conduction. In our simulation, we include radia- 
tive cooling, using t he st andard cooling function from 
ISutherland fc Dopital (| 19931 ) , conveniently extended towards 
higher temperatures, T > 10 7 K, which are attained within 
the cocoon generated in this simulation (see App. A in Paper 
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Figure 1. Temporal evolution of the star formation rate, for different snapshots of the simulation. In the upper panel we show the 
SF as a function of time (measured in Mq /yr and Myr respectively), marking with numbers 1, 2, 3, 4 the SFs at the times t PS 
0.4, 1.0, 1.2, 2.5 Myr respectively. In the four bottom panels we show the propagation of the cocoon within the ISM and the effect on 
clouds (corresponding to the numbers in the top panel). SF rate is obtained using the Schmidt-Kennicutt law (ra = 1.4) and collecting 
the clouds within a radius of 27 kpc. At very late epochs (t ~ 20 Myrs) the SFR declines to a very small value, as seen in the top panel. 
Note that our simulation maps only one-half of the galaxy, thus, assuming symmetry, we multiply our SF by a factor 2. 
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I). We also take into account gravity, but we neglect thermal 
conduction and large-scale, ordered magnetic field^]. 

The jet is modelled as a one-component fluid, charac- 
terized by a density pj which is a fraction e,- of the initial 
density of the interstellar medium. In order to prevent nu- 
merical instabilities at the jet/ISM injection interface, we 
use a steep, continuous and differentiable tra nsverse velocity 
and d ensity profile, as we did before (paper I; iPerucho et al.l 
12004 2005): 



V, 



cosh 



{(*?)"} 

cosh {(^)"} 



(i) 



(2) 



where: oij = 10 is an exponent which determines the steep- 
ness of the injection profile, rij, n env denote the jet and en- 
vironment electron number densities and the scalelength dj 
characterizes the width of the jet. The power injected by the 
jet is then given by Pj = P2nd 2 jPjVf, with f3 ~ 0.7158. In 
this simulation we set: Pj = 10 46 ergs -1 . 

We model the environment, where the jet propagates, 
as a 2-phase interstellar medium, comprising a hot, dif- 
fuse, low-density component and a cold, clumped system 
of high density clouds in pressure equilibrium with the dif- 
fuse component. The warm phase is characterized by a den- 
sity profile Penv(r) and a constant temperature T env . We 
assume that the diffuse gas is embedded within a dark 
matter halo, the latter being described by a NFW den- 
sity profile. This DM halo is chosen to have a total mass 
Mh = 5 x 10 11 Mq, concentration c = 10.2 and a scale- 
length l h = 206 h' 1 kpc. The ISM gas is assumed to be a 
fraction m g = Q ga s /CIdm ~ 0.212 of the DM, and to be 
distributed in hydrostatic equilibrium within this D M halo, 
follow ing the prescription given in Appendix C of iHesterl 
(2006). Each of the clouds in the cold comp onent is mod- 
elled as a truncated isothe rmal sphere (TIS: IShapiro et al.l 
Il999l . llliev fc Shapiroll200ll ). because this model seems to ad- 
equately reproduce the properties of clouds formed in sim- 
ulations of a thermally unstable ISM. TIS spheres posses 
a finite radius, and are characterized by two parameters, 
which we assume to be the mass M c i and a typical radius 
r c i. They are exact solutions of the equilibrium equations 
for isothermal spheres confined by an external pressure. We 
distribute 300 clouds within the simulation volume, using 
a mass spectrum previou sly derived from numerical simu- 
lations l|Baek et al.ll2005T ). Summarizing, the simulation is 
specified by 10 parameters: three of these describe the DM 
halo (Mh, c, lh), two the diffuse phase of the ISM (m g ,T env ), 
two the mass distribution within the cold component of the 
ISM (M c i,r c {), and the last three describe the jet (rij ,dj , Pj). 

We choose a simulation box having a size: Ltox = 
40 h" 1 kpc, so that the jet will diffuse through it at the end 
of the simulation. The spatial resolution attained is a func- 
tion of the maximum refinement level and of the structure 
of the code. For a block-structured AMR code like FLASH, 



1 Even an initially weak, small-scale, tangled magnetic field 
would be amplified to a level which inhibits thermal conduction, 
by 2-3 orders of magnitude w.r.t. the classical Spitzer values. 



Table 1. Parameters of the simulation. The halo mass is in units 
of Mq, distances are measured in kpc h , r ^env,^j are, respec- 
tively, the central gas electron and jet's densities, in units of cm -3 , 
T env is in Kelvin, dj and Pj are the jet width and power, respec- 
tively. The latter is expressed in cgs units (ergs- sec. 
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Figure 2. Star formation rate in 14 circular, equally-spaced 
shells. The results for different times are shown, from continu- 
ous to long dashed lines spanning t = 0.4, 1, 1.6, 2.5, 3.2 Myr. 



mum resolution along each direction is given by L^ ox /(n x 2 l ), 
where I is the maximum refinement level. In this simulation: 
= ri y = 8 and I — 6, thus the minimum resolved scale 
is 78.125 pc. Note that we are performing a 2D simulation, 
but we do not impose any special symmetry. 

We assume that the duty cycle of the jet is active 
for a time of 4 x 10 6 yrs , as i s typical for jets having 
such power l)Shabala et al.l [2008) . and its power declines 
linearly with time since this epoch until it switches off at 
t ~ 2 x 10 7 yrs. As for the star formation rate, we as- 
sume that the clouds are converting gas into star s with a 
rate specifie d by the Schmidt -Kennicutt (SK) law (jSchmidtl 
1 19591 . 1963, iKennicutd 119981 ): SFR = ± = AT, n , where 
A = (2.5 ± 0.17) x 10" 4 M Q yr _1 kpc -2 and n = 1.4 ± 0.15. 
We further require that SF is ongoing only within r egions 
having a mass larg er than the Bonnor-Ebert mass (Ebcrt 
1 19551 . lBonnor| [l956'). and temperature less than a specified 
upper limit, i.e. T ^ 1.2 X 10 4 K. These are rather conserva- 
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tive constraints, which tend to underestimate the magnitude 
of the negative feedback. 

Our initial setu p is very different from that o f the very 
recent simulation bv lSutherland fc Bick ncll (200?)), particu- 
larly in one very important point: while we model a popula- 
tion of pressure-confin e d clou ds embedded within the ISM, 
ISutherland fc Bicknelll ([2007) put a turbulent disk around 
the jet's source, embedding the AGN within it, as they are 
more interested in addressing problems related to the evo- 
lution of GPS and CSS radio sources. Also, note that their 
simulation box is much smaller than ours (1 kpc), and the 
temporal scale is also very different (« 10 5 yrs in their sim- 
ulation, compared to 2 x 10 7 yrs in the present work). More- 
over, one may notice that we have neglected the velocity 
field of t he inhomogeneous component : in fact, as already 
noted bv lSutherland fc Bicknenl (|2007h . for the spatial and 
temporal scales of interest, its inclusion has little effect on 
the evolution of the clouds and of their SF. The typical tur- 
bulent velocities of clouds are in the range 10 — 70kms _1 . 
Thus, on typical timescales <^ 2 x 10 7 yrs in our simulation, 
these clouds have displacements I < 0.2 — 1.4 kpc, smaller 
than the average distance among clouds (see Fig. [1} ■ 



3 PROPAGATION OF THE JET 

Soon after the jet enters the ISM, a low-density region, 
the cocoon, is generated. The details of the propagation 
of the jet within the ISM, and the properties of the co- 
coon , have been extensively studied in the past |Scheuer] 
1 19741 . Fallc 199l]). Only recently, however, numerical sim- 
ulations have been used to study in more detail the effect 
which this interaction has on the inhomo geneous component 
of the ISM, where stars are forming fe.g lSaxton et al1l2005l . 
iKrause fc Alexa nder 2007|). In Paper I, we extensively ana- 
lyzed properties of cocoon and the interaction with a single 
cloud that is forming stars. In the following we will make few 
comments about cocoon propagation, concentrating rather 
on the effects on our multicloud system. 

3.1 Quenching of star formation 

The global evolution can be seen in Fig. [T] it is possible to 
distinguish two main phases, corresponding to the evolution 
of the cocoon: an active and passive phase. While the jet is 
active, the cocoon is fed and it expands almost self-similarly, 
while, when it is switched off, the cocoon diffuses and affects 
a larger fraction of the ISM. 

Initially (see snapshots 1 and 2 in Fig. [TJ, the inter- 
action mostly takes place through the weak shock at the 
interface between the cocoon and the ISM. This shock com- 
presses and heats up the clouds, having two opposite effects 
on SF: compression tends to increase SF, but the larger tem- 
perature also increases the critical mass for gravitational 
collapse, thus decreasing the volume fraction of the clouds 
which can actively form stars. Globally, we observe an ini- 
tial enhancement of SF (i.e., we observe a positive feedback). 
The effects of the two concomitant processes on clouds can 
be seen in Fig. [2] where we build the shell-SFs by adding the 
SF of clouds lying in a concentric annulus of mean radius 
r and plot them as a function of radius for different time 
steps. Here, few episodes of positive feedback are evident 
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Figure 3. Timescale for feedback. The three curves refer to the 
same density profile, having po = 1-7 X 10 -25 gcm — 3 , /3 = 2, 
and 3 values for ao = 0.8, 1, 1.2 kpch (dashed, continuous 
and dotted curves, respectively). The strong dependence on ao is 
simply a consequence of the high exponent with which it appears 
in the expression for tQ. 

in the clouds which are nearer to the cocoon. Later ( see 
snapshots 3 and 4 in Fig. [TJ , the cocoon propagates within 
the medium inducing a general increase of the temperature, 
which heats up the clouds and decreases their density, and 
finally the outer regions of the clouds are stripped due to 
Kelvin-Helmholtz instabilities, as we studied in more detail 
in Paper I. These effects tend to reduce the mass of the 
clouds, decreasing SF, and eventually drastically suppress- 
ing it on a time-scale of 2 — 3 x 10 6 yr (negative feedback). 
Obviously, the more external regions in the galaxy are af- 
fected later by the disruptive effect of the cocoon, being 
almost unperturbed until t ~ 10 6 yr and t = 1.6 x 10 6 yr, 
respectively at r ~ 30 and ~ 45 kpc (see Fig. [5J. As we 
see from this simulation, approximately on a timescale over 
which the cocoon diffuses, the clouds are destroyed: this is 
then the typical time-scale over which SF will be inhibited. 
Obviously, one must observe that this timescale is not ex- 
actly coincident with the duty cycle of the jet. 

Following Paper I, we restricted our simulation to a sin- 
gle episode of jet injection. Any subsequent event of jet in- 
jection would have little influence on SF, since the jet would 
propagate through a high temperature (T ~ 10 8 — 10 11 K) 
and low density (n e ~ 10~ 2 — I P " 1 cm " 3 ) environment. As 
already noted in llnoue fc Sasaki (|200ll ) and in Paper I, after 
the injection of the jet, the cooling time of the diffuse ISM 
exceeds the dynamical time by a factor 6.5 x 10 2 — 3 x 10 5 ; 
thus the heated gas is not able to cool, quenching the second 
emitted jet. The extent of this region crucially depends on 
the physical parameters of the medium and on Pj . The case 
we present in this paper is that of a very powerful jet, so at 
the end all the gas within the affected region is influenced 
by the expansion of the cocoon. For less powerful jets, it 
is reasonable to expect that the extent of the quenched SF 
region will be much smaller. 

We should however remark that this scenario could not 
be exhaustive of all the possibilities. In a recent paper, 
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Figure 4. Estimated synthetic colours FUV-NUV, NUV-g and g-r (in AB system) as a function of galaxy age (in units of Gyr) for 
Z = Zq. Red, blue, green and orange curves correspond respectively to an unperturbed exponential SF with r = 0, an exponential SF 
with t = 1 Gyr and *AGiV = 1 Gyr, an exponential SF with r = 1 Gyr and tytGiV = 5 Gyr and an unperturbed exponential SF with 
r = 1 Gyr. In the inset panels we show the difference of colours with respect to the reference one with an unperturbed exponential SF 
with t = 1 Gyr as a function of galaxy age; the colour code is the same as in the plot of colours vs age. 



iKrause fc Alexander! (|2007l ) studied the evolution of clouds 
lying very near to a jet. The clouds are destroyed by the 
Kelvin-Helmoltz instability produced by the jet, but the au- 
thors notice that a fraction of these clouds, under the effect 
of thermal instabilities, evolve into filamentary structures 
which are characterized by high densities. These filaments 
could then still host SF, although it is not clear the quanti- 
tative relevance of this phenomenon to global SF within the 

host g alaxy. 

In lSutherland fc Bicknelll (|2007h the mechanical energy 
of the jet is mostly dissipated when the jet tries to perco- 
late through the disk, while when the jet has already formed 
a cocoon, as in our simulation, its energy feeds the turbu- 
lence within the cocoon. On the large scales probed by our 
simulation, the turbulence within the cocoon is responsi- 
ble for the final destruction of the clouds. When the jet is 
switched off this turbule n ce wi ll decay but, as already no- 
ticed by llnoue fc Sasakil (|200lh . often the cooling time of 
the very hot plasma is too long to allow the formation of 
cold clouds. Other external episodes like minor interactions 
between clouds could however fuel again some cold gas, and 
strip the hot gas, thus favoring the subsequent development 
of embedded cold clouds. 



3.2 Time-scales 

Although we discuss the output of only one single simula- 
tion, the general model of jet propagation in the ISM which 
is also probed by our simulation allows us to determine one 
of the most relevant parameters of the star feedback model 
that we will develop in the next paragraphs: the typical time- 
scale for suppression of stellar formation. 

In the se lf- simila r expansion model of 
iKaiser fc Alexander! (|l997T I. the cocoon is supposed to 
propagate into an unperturbed ISM which is well described 
by a power-law profile: p(r) = po(r /ao) _/3 , and the typical 
scale length of the jet varies with time according to: 

/ t \ 3/5-/3 

L J =c 1 a ^—J (3) 
where (|Kaiser fc Alexander! [l997l , Eq. 5): 



to = 1.186 x 10 6 (j|^y /3 yrs. (4) 

and Pj,45 is the jet's mechanical power in units of 
10 45 ergs" 1 . 

We have assumed that the number density of cold, star- 
forming clouds is proportional to that of the diffuse gas, and 
thus can be well approximated by the same power-law den- 
sity profile outside the core (r > an). There exists a thresh- 
old cloud number density under which SF becomes negligi- 
ble, and we will suppose that this corresponds to a value of 
the gas density p cr « 10 -27 gcm~ 3 , which is reached at a 
distance: r cr = ao(pn/p cr ) 1/ '' 3 ■ Inserting this into Eq. ©, we 
find that the characteristic timescale needed to reach this 
distance is given by: 

t/6 = to f c x (5) 

We plot in Fig.[3]the dependence of tfi, on jet's power Pj. As 
we can see, typical values for the time-scale of feedback are 
in the range 10 6 - 10 s yrs, for reasonable values of the input 
parameters. The values of Pj that we choose are represen- 
tative of the observed mass range of BHs at the centers of 
typical galaxies, accord ing to the empirical relat ion between 
M bh and Pj found by iLiu, Jiang, fc Gul (120061 ). Thus, the 
quenching time-scale for our simulation, of the order of a 
few Myrs, will increase of about two orders of magnitudes 
for less powerful jets. 



4 SYNTHETIC COLOURS 

The simulation discussed in this paper describes the impact 
that the jet emitted by AGN has on quenching SF. However, 
we would like to transpose this modification of the SF into 
a change of some observable quantities. Colors and absorp- 
tion/emission lines are the most important observables given 
by the observations: here we will concentrate on the former 
ones. In order t o reproduce galaxy c olour s, we use the syn- 
thetic models of lBruzual fc CharlotJ (|2003l . hereafter BC03), 
that allows us to encompass a wide range of metallicities, 
starting from galaxy ages t ga i of 10 5 yrs, and gives a full 
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Figure 5. Estimated synthetic colours FUV-NUV, NUV-g and g-r (in AB system) as a function of galaxy age (in units of Gyr) for 
Z = Zq and for r = 1 Gyr (red lines) and r = 3 Gyr (blue lines). Continuous, long-dashed and short-dashed lines correspond to 
t AG N = 1, 5, 10 Gyr. 
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Figure 6. Estimated synthetic colours FUV-NUV, NUV-g and g-r (in AB system) as a function of galaxy age (in units of Gyr) for 
r = 1 Gyr and tAGN = 5 Gyr. Long-dashed, continuous and short-dashed lines correspond respectively to Z = 0.008, 0.02, 0.05. 



coverage in wavelength from 91 A to 160 /im. We start from 
a single initial burst model, consisting of a single population 
with a Chabrier IMF0 and three values of the metallicities: 
Z = 0.008, 0.02 i.e. (Z = Z e ), 0.05 and convolve it with a 
SFR law. 

We suppose that without interaction of the AGN with 
the ISM, SF within the clouds evolves following an expo- 
nential SF law SFR oc e~*' T , where r is a characteristic 
timescale corresponding to the time when SF is reduced of 
e .. Then at t = tAGN, the AGN interaction begins to 
work, shocking the ISM and inhibiting SF. This configura- 
tion can be realized by combining the unperturbed expo- 
nential SFR for t ^ tAGN and SFR from our simulation 
for t > tAGN- Due to the short time-scale involved in the 
AGN effect on SFR, positive feedback is not observable in 
the derived colours, while the negative feedback is directly 
translated into a strong truncation of SFR. Note that this 
model is an approximation, since less powerful AGNs are 
observed at low redshifts and the decrease of SF can be 
shallower, as discussed in the previous section (see, also, the 
parame tric model of SF quenching discussed in lMartin et al.l 
|2007al )). If t — » co, the SFR becomes a burst with finite 

2 Galaxy colours are unchanged if we use a Salpeter IMF, since 
these two IMFs differently describe the distribution of low mass 
stars that contribute little to th e light distribution, w hile strongly 
affecting the total stellar mass dTortora et alj|2009ft . 



length tAGN (i.e. a SF constant up to tAGN with a null value 

for t > tAGN)- 

We calculate colours by convolving the filter responses 
of u, g, r, i and z SP SS and the ultrav i olet FUV and 
NUV GALEX band s (|Martin et all 120051 . lYi et all 120051 , 
iKavirai et al.l [2007a) with synthetic spectra. In Figs. |31 O 
and[6]we show the change with time of colours FUV-NUV , 
NUV — g and g — r and the effect of quenching. In partic- 
ular, in Fig. [3] we compare the colours of our inhibited SF 
(with t ag N = 1 and 5 Gyr) and two simple models using an 
exponential SF with r = and r = 1 Gyr. SF inhibition es- 
sentially reddens the colours starting at t>tAGN- Ultraviolet 
is much sensitive to SF inhibition (see the middle panel in 
Fig. [4| , with variations of 1 — 3 magnitudes in NUV — g; less 
sensitive are the visible colours, with variations of ~ 0.2 — 0.4 
in g — r. The colour FUV-NUV increases after the action of 
the AGN up to ~ 1 — 4, depending on tAGN, but after this 
event it decreases reaching a value < 1 for large t ga i . There- 
fore, the latter colour is degenerate, because one is not able 
to univocally determine the galaxy age from colour (each 
colour corresponds to two possible values of t ga i). 

A single burst model (r = 0) is similar to our models 
with inhibited SF for t > 2 — 4 Gyr, but predicts larger values 
for t < 2 — 4 Gyr. While g — r colours for the different models 
shown in Fig. U are similar since t ~ 6 — 8 Gyr, the other 
colours are comparable at larger ages t ~ 8 — 10 Gyr, being 
more sensitive to SF. Not shown in the figure is the case 
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Figure 7. UV/optical colour-colour diagram (NUV-r vs g-r) from synthetic spectra redshifted to median redshift of the sample z me( ; = 
0.057 assuming r = 1 Gyr and metallicitics Z = 0.008 (long-dashed line), Z = 0.02 (continuous line) and Z = 0.05 (short-dashed line). 
Red, blue, green and yellow lines correspond, respectively, to <AGiV = 0.1, 1, 5, 10 Gyr. The large points set the values of t^GN on 
each synthetic track, while the small ones indicate t ga i (= 0.1, 1, 5, 10 Gyr), which on different tracks correspond to different set of 
colours after the quenching. In the inset panel we superimpose the sample galaxies to theoretical tracks, with violet, gray and cyan 
points indicating galaxies in the redshift bins 0—0.04, 0.04-0.06 and 0.06—0.08. Instead of error bars, we show in the right panels the 
distributions of uncertainties on galaxy colours g — r and NUV — r. 



t — > oo, characterized for fixed values of tAGN and t ga i by 
bluer colours. 

These considerations are obviously dependent on the 
choice of r and Z. In Fig. [S] we compare the results for 
t = 1 and 3 Gyr. At fixed age, a more protracted un- 
perturbed SF predicts smaller NUV — g and g — r and a 
larger FUV — NUV. Metallicity has the opposite effect: 
FUV — NUV decreases if we increase metallicity, while g-r 
increases; finally, a more complex behaviour is observed for 
NUV — g, as can be seen by inspecting Fig. [U These com- 
plex dependencies on metallicity are linked to the details of 
stellar population prescriptions. 

We can outline a picture where the AGN effect has a 
main role in the evolution of brighter elliptical galaxies. In 
particular, inhibition of SF by AGN transforms a galaxy 
with a protracted SF (in our case an unperturbed SF ob- 
tained using an exponential SF law with a specific r) into a 
more quiescent galaxy with SF stopped at a time approxi- 
mately equal to the epoch of jet injection. The extreme case 
of this model is represented by a single burst model, where 
only at t = is the SF observable and most of the SF has 
been completed by this epoch. 



5 EPOCH OF QUENCHING EVENT 

We will now attempt to fit synthetic models to a large sam- 
ple of local ETGs, in order to obtain information about the 
epoch of quenching event and RSF. The galaxy sample is 
presented in Sect. [570 while the spectral fit procedure and 



the first results are described in Sect. 15.21 We go into more 
detail in Sects. [5^3l and [5^4l where quenching and the proper- 
ties of the recemtly formed stellar populations are discussed. 
Finally, few comments on galaxy evolution are addressed in 
Sect. [53] 



5.1 Galaxy sample 

We use a sample of ETGs extracte d from the SPSS, us - 
ing a selection procedure described in lKavirai et all (1200731 ). 
The initial selection is made using the fracDev parameter in 
SDSS, that attributes a weight to the best composite (de- 
Vaucouleur's + exponential) fit to the galaxy image in a par- 
ticular band. The criterion fracDev > 0.95 has been proven 
to be extremely robust, allowing one to pick up ~ 90% of 
ETGs in a typical sample of SDSS galaxies. Furthermore, 
a visual inspection of SDSS images is needed to refine the 
selection: the ability to classify galaxies obviously depends 
on redshift and apparent magnitude of the observed galax- 
ies. In order to construct a magnitude-limited sample, we 
restrict ourselves to r-band magnitude < 16.8 and redshift 
2 < 0.08. Finally, cross-matching with ultraviolet GALEX 
data produces the final sample, with measured magnitudes 
in S PSS bands u g, r, i z and GALEX FUV and NUV 
(see iKavirai et all l|2007al ) for further details). In this way, 
we have a wide coverage of galaxy spectra, since in the opti- 
mal cases, magnitudes in 7 bands are observed, ranging from 
A ~ 1500 A up to A ~ 9000 A. SPSS magnitudes have typical 
uncertainties of ~ 0.01, while GALEX data are more uncer- 
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Figure 8. Distribution of recovered best fit values of t ga i — tAGN 
for different spectra sub-libraries extracted from our initial li- 
brary. On the y-axis, we plot the fraction of galaxies (normalized 
to dimension of the sample) within bins of t ga i — tAGN- I n the 
range t ga [ — tAGN 6 [0 — 2] Gyr the scale is linear, while for 
tgai ~ tAGN > 2 Gyr we use an arbitrary logarithmic scale, fn 
the linear range we group galaxies in bins of size 0.2 Gyr, while 
larger bins are used outside this range. In particular, we show 
the results leaving all parameters free to change, setting Z = Zq 
and T = I, 3 or 5 Gyr, for r = 3 Gyr we analyze two cases with 
constraints on tAGN and then include internal extinction in the 
fit. We analyze the case of a finite length burst (r = oo) and 
a combination of two coeval stellar populations with r = and 
t — 3 Gyr with a free mass ratio or a fixed one. See labels in the 
plot and discussion in the text for further details. 



tain with mean errors of ~ 0.25 and ~ 0.15, respectively for 
FUV and NUV magnitudes. 

In Fig. [7]we show a colour-colour (NUV — r vs g — r) dia- 
gram from our synthetic tracks. In the inner plot we present 
our sample, superimposing the galaxies, colour coded ac- 
cording to redshift bins. This diagram gives information 
about both t ga i and tAGN selecting specific values for two 
colours and the best parameters of galaxies. Different values 
of galaxy parameters predict galaxy colours that populate 
different regions of the diagram. 

The sample of lKavirai et alj (|2007aT ) is suitable for our 
analysis, since AGN spectral features are observed in spectra 
of many galaxies. Type-I AGN are automatically removed 
by using the SDSS spectral classification algorithm, while 
we are interested in galaxies that host a Type-II AGN. To 
distinguish between normal star forming galaxies and galax- 
ies hosting an AGN, it is usual to analyze a few strong 
emission lin es, e.g. the emission li ne ratios \OIII /H 8] an d 
[Nil /Ha] (|Baldwin et all Il98ll . iKauffmann et al. 1 120031 ). 
Galaxies with both these indicators measured represent 
~ 65% of the galaxy sample, and ~ 86% of them have 
spectral features consistent with those of LINER, Seyfert 
or transition objects. 

Emission from Type-II AGN does not affect the stel- 
lar continuum of host galaxies. In fact, for luminous Type- 
II AGNs the maximum contamination in flux amounts to 
few percents in visible bands and to less than 15% (trans- 
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Figure 9. Median K-corrections derived by fitted synthetic spec- 
tra for bands FUV, NUV, r and colour FUV-NUV and NUV-r. 
The bars are the median scatters. See legend for details on the 
colour code. 

lating to 0.15 ma gs) in UV bands (|Kauffmann et al~ll2003l . 
ISalim et al.ll2007l) . In addition, galaxies hosting AGN are 
systematically redder in the UV colours than their coun- 
terparts which do not have AGNs, further on suggesting 
that continuum emission from AGN leaves unaffected galaxy 
spectra. 

5.2 Spectral fit 

We build a library of synthetic spectra, using our SFR pre- 
scription with galaxy age t ga i, r, tAGN and Z as free pa- 
rameters. We divide the sample into 3 redshift bins 0-0.04, 
0.04-0.06 and 0.06-0.08 and we move our spectra to the 
median redshift of each bin (respectively of 0.031, 0.052, 
0.069). Synthetic magnitudes and colours are obtained con- 
volving these redshifted spectra with the filter responses. 
Finally, the synthetic colours are fitted to the observed ones 
(FUV - NUV, NUV -g,u-g,g-r,g-ia,ndg-z),by 
a maximum likelihood method, which allows us to estimate 
the best values for the free galaxy parameters. In this way 
we do not need to have K-corrections, leaving the fit safe 
from possible uncertainties introduced by these corrections, 
and the simplistic division of our sample into only 3 redshift 
bins will not affect our estimates. Later, we will discuss the 
K-correction that we derive from the fitting procedure and 
we will use to obtain rest-frame magnitudes. 

We restrict synthetic library to spectra with t ga i free to 
change (with a tiny step) up to 17 Gyr, tAGN £ (0, 15) Gyr, 
Z e (0.008, 0.02, 0.05) and 3 reference values for the SFR 
scale t = 1, 3, 5 Gyr. For each value of t ga i we have more 
than 150 spectral models. This libr ary is wide enough to re- 
produce spectral features of ETGs l|Panter et al.ll200&t ). The 
range of metallicities used here has been show n to be repre- 
sentative of luminous ET Gs with Mb < — 19 (|Romeo et alj 
120081 . iTortora et ail 2009). In particular, the fit of spectra 
with an unperturbed exponential SF and variable Z to lo- 
cal ETGs gives on average values of r < 1 and Z > Zq 
l|Tortora et all [2009 ) . Such an unperturbed exponential SF 
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Figure 10. Colour-magnitude diagrams of our sample, colour-coded following the classification shown in the legend: t gat — t AG N < 
(orange), < t ga i - t AGN < 0.5 (green), 0.5 < t gal - t AGN < 1 (blue), 1 < t gal - t AGN < 4 (violet), 4 < t gal - t AGN (red). We 
show 2 = colours predicted using a Zf = 3 single burst population for Z = Zq as a continuous gray line and various Z < Zq (i.e. 
Z = 0.0001, 0.0004, 0.004, 0.008) as dashed gray lines. Left Panels. We plot g-r vs r and u-r vs r diagrams. Right Panel. We plot NUV-r 
vs r magnitude. 



with a low time-scale r can disguise a more complex SF 
evolution, which we have obtained from the simulation and 
want to probe against observations. 

Before proceeding to model the properties of the ob- 
served sample of galaxies, we briefly discuss the possible sys- 
tematics which the fitting procedure can generate in the best 
fitted parameters, by performing a set of Montecarlo simula- 
tions on synthetic colours. We have extracted a large sample 
of simulated spectra from our SED library with random t ga i, 
tAGN, t and Z. Then, we applied our fitting procedure and 
compared the recovered best fit parameters against the in- 
trinsic ones. While t ga i, tAGiv0 and Z are recovered quite 
well, r is poorly constrained. To reduce the unavoidable de- 
generacies in the fitting procedure, it would be acceptable 
to set r to a fixed value, however, in the following we will 
discuss both the cases with r variable and constant. In addi- 
tion, to further reduce the noise in our final results we also 
analyze the effect of adopting constant metallicity. 

Coming back to the observed sample of galaxies, we 
use different synthetic libraries to fit the observed colours, 
extracting spectra from our initial and extended set of syn- 
thetic models. We show the results of this analysis in Fig. 
[8] where we plot the distribution of recovered values of 
tgai — t AG N . If on one hand the change of spectral library can 
modify the estimates of single values of the parameters, on 
the other hand it leaves the main clump of the distribution 
of differences t ga i — t AG N unaffected. Due to the correla- 
tions among parameters and constraints imposed on some of 
them, the sample distributions of recovered values for t ga i, 
t AG N and Z can change, but median values for t ga i — t AG N 
(~ 0.8 - 0.9 Gyr with a median scatter of ~ 0.1 - 0.2 Gyr) 

3 Note that the input value for t A GN is correctly recovered if 
tgal — tAGN 0, while for t ga i — t AG jy < it is not possible to 
constrain its value. 



are almost unchanged, with a scatter among the different 
estimates consistent with 0. If we impose strong constraints 
on t AG N (e.g., t AG N ^ 5 Gyr), we observe a departure from 
the mean trend obtained using other libraries, with a peak at 
very large values of t ga i —t AG N (corresponding to quenching 
at high redshift), but this constraint does not seem to be mo- 
tivated since it decreases the quality of fit. In addition, note 
that for t = 1 Gyr the distribution peak is slightly lower and 
less prominent, mainly due to the larger number of galaxies 
estimated to have unperturbed SF (i.e., with t ga i < t AG ^). 
To perform a more complete analysis, we also show the re- 
sults for a combined spectral model obtained by summing 
up two coeval stellar populations: a single burst model to 
a quenched SF assuming a variable or a fixed mass ratio of 
the two populations. This model takes into account the co- 
existence of stellar populations having different properties: 
a part of stars formed in a single initial burst, in addition 
to this component, there is gas which cools to form continu- 
ously stars and is affected by AGN. What we see in Fig. [8] is 
that these combined models we have discussed give the same 
results as using the single component ones. Our estimates of 
tgai — t AG M are robust, since the bulk of their distribution 
depends little on the duration of the unperturbed SF and 
Z. Despite this result, as we will see later, the past SF his- 
tory of galaxies is dependent on the choice of these models. 
Finally, two shortcomings have to be discussed. If consider- 
ing more protracted unperturbed SF we will obtain a larger 
number of galaxies with higher t ga i — t AG M values than using 
SFs with lower r (see Fig. [Sj . In addition, a similar effect 
is given if quenching is not almost instantaneous as in our 
simulation, also in this case we would obtain larger values 
of tgai — tAGN- However, our fast quenching is a good ap- 
proximation to describe AGN effect for a wide sample of 
galaxies, but further analysis on these aspects waits to be 
done in future. 
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Figure 11. UV/optical colour-colour diagram. We plot NUV-r 
vs g-r of galaxies in our sample colour-coded following the classi- 
fication shown in the legend. 



In order to avoid confusion, we will use in the following 
the results obtained in the most general case where Z, r, t ga i 
and tAGN are all unconstrained. Also, K-corrections deserve 
particular attention. In fact, due to uncertainties in the UV 
spectrum, it can be difficult to quantify the amount of this 
correction, acircumstance which can introduce systematics. 
In our range of redshifts, we find typical NUV corrections of 
— 0.3. in good agreement with those found in lKavirai et al.l 
l|2007al ) using the same galaxy catalog. In Fig. [9j we show 
also results for K-corrections of FUV, r and colours FUV- 
NUV and NUV-r. We obtain a sharp rise in NUV correction: 
this is negligible in the lowest redshift bin, while its median 
value becomes ~ 0.25 at higher z. The median scatter in the 
estimated values of NUV and FUV corrections is consistent 
with typical unce rtainties of ob s erved m agnitudes in these 
bands. Note that iKavirai et al. (12007a ). using a 9 Gyr old 
single burst model and lRawle et al.l ( 20081 ). using both single 
burst and frosting models, find larger corrections of 0.2 — 1.2 
for galaxies having z <^ 0.1. 



5.3 Quenching event 

In Fig. [TU]we show the visible and UV colour- magnitude di- 
agrams for our sample. As already detailed in lKavirai et al.l 
(2007a), the UV dia gram in the right panel of this figure 
is largely different from similar colour diagrams for optical 
bands (see left panels), since it shows a larger scatter. Fig. 
1101 shows the power of NUV — r colour in selecting galax- 
ies with different SFR histories. In particular, the points in 
the figure are colour-coded according to the fitted value of 
tgai — tAGN'- galaxies that experienced some quenching phe- 
nomena in the past lie upper in this diagram, while those 
having lower values of this quantity populate on average a 
region with a flatter UV continuum and a larger scatter. 

A rough comparison with the expectation of the mono- 
lithic scenario has been made, superimposing the predicted 
colours for single burst populations with a more extended 



Figure 12. t ga i — tAGN as a function of ultraviolet magnitude 
NUV, while in the inserted panel is plotted t ga \ — tAGN as a 
function of colour NUV-r. Points are colour coded for various 
redshift bins. See Fig. \E\ for details on the scale used in the y- 
axes. 



formation redshift starting at zf — 3. For completeness, in 
addition to a model with Z = Zq we also show the pre- 
dictions for various metallicities < Zq implemented within 
BC03 prescription. These predictions are not very sensitive 
to the precise value of z; (at least when considering rela- 
tively old stellar populations), and our choice is consistent 
with optical analyses which estimate Zf ;> 2 (jBower et al.l 
1992). We see that an almost solar single burst model is 
only consistent with redder galaxies both for g-r and NUV-r 
colours. Thus, while galaxies with low g-r and NUV-r can 
be matched by models with low Z, galaxies in the u-r vs 
r plot appear to depart from these predictions which sys- 
tematically lead to redder colours. However, as we pointed 
out in a previous subsection, luminous ETGs are almost 
solar and very low values of metallicities lik e those used 
in Fig. | 10 | are not applic able to our galaxies (jPanter et al.l 
120081 . (Tortora et all 12009). Moreover, semi- analytic simula- 
tions predict that a single so lar metallicity is a reasonable 
hypothesis to describe ETGs (Nagashima fc Okamotdl2006l ) 
and similar results for the brighter and more massive galax- 
ies have been obtained in N-body+hydrodynamical simu- 
lations (|Romeo et al.l 120081 ) . Thus, the comparison with a 
solar metallicity single burst model indicates that only a 
few galaxies (< 6 — 7%) are consistent (within the errors) 
with a single and old stellar population; these galaxies are 
the redder ones mainly with t ga l — tAGN }> 1. This result is 
co nsistent with the < 1% of purely passive galaxies found 
in IKavirai et al.l (|2008l ) analyzing a sample of galaxies at 
z J> 0.5, since in the local universe we observe a larger num- 
ber of quiescent and passive galaxies. 

In Fig. [TT] we show NUV-r vs g-r, with points colour- 
coded according to Fig. 1101 A correlation is observed, but 
note that the excursion in NUV-r is much larger than in g-r. 

Galaxies with a stronger ultraviolet flux are on aver- 
age characterized by higher RSF, that translates into lower 
values of t ga i — tAGN, while those with a lower NUV flux 
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Figure 13. NUV (and NUV-r) vs. RSF for a model leaving free all parameters (left panel) and one with Z = Zq and r = 3Gyr (right 
panel). The colour code is the same as in Figs. ITuland llll 




Figure 14. NUV (and NUV-r) vs RSF for a model leaving free all parameters. In the left panel we colour for Z, with Z = 0.008 (green), 
Z = 0.02 (blue) and Z = 0.05 (red). In the right panel we colour for t with r = 1 Gyr (red), r = 3Gyr (blue) and r = 5Gyr (green). 



have SF that is quenched early. This correlation is shown in 
Fig. 1121 In the inner panel, we show the trend with colour 
NUV-r, implicit in Fig. [TO] 



The results concerning the scale of SF are also inter- 
esting. When considering the results obtained fitting our 
reference spectral library, we find that only ~ 15% of galax- 
ies have r = 1 Gyr, while a more protracted SF is recov- 
ered for the other galaxies in the sample, with ~ 43% and 
~ 42% of these having r = 3 and 5 Gyr. As a compari- 
son, we also fitted unperturbed SFs with r 6 (0.1, 5) Gyr 
and Z g (0.008, 0.02, 0.05) obtaining that ~ 86% of galax- 
ies have r ^ 1, whi le less than 1% have r > 3 Gyr (see also 
iTortora et al.|[2009l ). Thus, the low SF scales recovered when 



simple unperturbed exponential SFa[3 are fitted to data can 
hide a galaxy population with a more protracted background 
SF that is quenched in the late stage of galaxy evolution for a 
feedback effect. These results are consistent with the general 
scenario depic t ed for the color evolution of E+A galaxies in 
iKavirai et all l|2007bl ). They find that, superimposed to an 
early burst of formation, a recent SF burst (which typically 
takes place within 1 Gyr) over a timescale ranging between 
0.01 and 0.2 Gyrs is needed to match galaxy colours. These 
galaxies are just migrating towards the red sequence and 
show a SF quenching which is correlated with their stellar 
mass and velocity dispersion and linked to different sources 



4 Similar considerations can be made if we assume a delayed SF. 
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of feedback (AGN or supernovae). However, although also 
different models can have the same final result as our un- 
perturbed and truncated SFs for different class of galaxies, 
this work reproduces typical timescales consistent with those 
discussed here (see also Fig. [3j . 

Fitting our reference library, we find that ~ 35% of 
galaxies have t ga i > tuniv, while this percentage decreases 
to ~ f 5% for those galaxies with t ga i > tuniv + 2. Excluding 
all galaxies with t ga i > tuniv we find a median formation 
redshift of Zf = l.Olto'4 (uncertainties are 25 th and 75 th 
percentiles); note that the distribution has a strong tail for 
high formation redshift. For tAGN we find the best estimate 
zagn = 0.13 ± 0.02. 

5.4 Stellar mass fraction and AGN feedback 

Until now, we have discussed how the quenching event is 
linked to colours and galaxy luminosities, but we have not 
expressly discussed the amount of stellar mass that is pro- 
duced. To quantify this recent star formation, we define the 
RSF via the amount of stars produced (i.e. the produced 
stellar mass fr action) in the last 1 Gyr in the rest frame 
of the galaxy dKavirai et al.|[2007al . 2008). As discussed in 
the previous subsections, this quantity is strictly dependent 
on the time of the quenching event, and in particular on 
tgai — tAGN- The latter has been show to depend on colour 
NUV-r and the luminosity in the NUV band, while it is less 
sensitive to the same quantities obtained only using visible 
bands. 

Galaxies with brighter ultraviolet fluxes and bluer 
NUV-r colours have much more RSF with respect to fainter 
ones. This result is not surprising, since an indication of it 
has been obtained by the trends shown in Fig. fTJ] In Fig. [T3] 
we show the RSF as a function of NUV and NUV-r. In the 
left panel, we show the results using our reference library. 
Obviously, only galaxies with t ga i — tAGN less than 1 Gyr 
survive in these plots and show the presence of RSF. In this 
case, the trends are not so tight, since some galaxies depart 
from them (red points). As a comparison, in the right panel 
we show the results for a model with fixed values of Z and 
r, where the correlations are tighter. Coming back to the 
results we obtained using our reference library, we show in 
Fig. Q3] the same results already plotted in the left panel of 
Fig. 1131 but now coded for metallicity and r, respectively in 
the left and right panels. These plots show that those galax- 
ies, which experience a significant RSF, also have metallic- 
ities significant different from the average. In fact, galaxies 
with different metallicities seem to stay on three different se- 
quences, with a lot of supersolar galaxies having red colours 
and a low NUV flux, but a large SF. This trend is more evi- 
dent if we look at the plot as a function of colour. The trends 
with SF time-scale r are less clear, but some of those galax- 
ies departing from the mean trend have r = 5 Gyr. Thus, 
some systematics might be introduced into these results by 
the degeneracies between Z and r, and possibly with other 
parameters. However, our results are robust, since qualita- 
tively we recover some trends which do not depend on the 
particular choice of spectral library. 

For the model that leaves all of the parameters free, we 
find a median RSF of ~ 0.2lg 3%. In particular, galaxies 
with the bluest colours (and with larger NUV fluxes) have 
RSF of ~ 1%, while the redder ones have RSF less than 



0.1%. On the contrary, a lower median RSF of ~ 0.1% is 
observed for the model with r = 3 Gyr and Z = Z@ . 

On average, our recove red RSF fractions ar e slightly 
lower than those obtained in lKavirai et alj (|2007aT ) . but this 
is not surprising since we follow a different approach. In de- 
tail, we adopt what we can call a blue-to-red approach, since 
we have modelled galaxies using an exponential SF. Thus, 
a galaxy is initially blue, and later it becomes red since SF 
is quenched. On the contrary, using a red-to-blue approach, 
galaxies are modelled with a single burst population, which 
predict systematically redder colours than those obtained 
with a more complex SF. Thus, to match the observations 
one needs a recent burst to make the colours bluer. It is how- 
ever worth to be noted that within scatter the two analysis 
give fully compatible results. 

Our results can also be connected with the downsizing 
scenario. In fact, the fact that redder galaxies (i.e. with a 
lower NUV flux) have less RSF with respect to t he bluer ones 
can reproduce the trend s found in recent wo r ks rtCowie et al 



2006, 



iBundv et al 



2006. 

-nr 2 



De Lucia et al 



1996. Borch et al 

200rj. iTragerl I2000L iThomas et al.l 120051. iNelan et al.| 1200 
Kavirai et aljboosl . [Romeo et alJl200St , iTortora et alj|20Q ~ 
However, this connection has to be better analyzed and fur- 
ther information could be recovered by extending the lumi- 
nosity range. 



5.5 SED evolution 

In the previous subsection, we obtained results that indicate 
the presence of a RSF that is connected with the quench- 
ing events. Here, we quantify the evolutive path of galax- 
ies in our sample, discussing how colours and magnitudes 
(averaged over the galaxy sample) evolve. A more detailed 
analysis of these results is beyond the scope of this paper. 

We verified that a large part of the galaxies in our sam- 
ple are not very sensitive to changes of the details of the 
spectral library, since they predict a median t ga i—tAGN that 
does not change siginficatly (see Fig. [8} . Thus, our models 
are able to successfully describe the recent formation his- 
tory of galaxies, while the extrapolation to early phases of 
galaxy evolution is influenced by the choice of unperturbed 
synthetic spectra. However, this extrapolation can fail to 
describe the earlier evolution, since galaxies which have ex- 
perienced various (minor or major) mergers, and AGNs ca n 
act to quench SF at different epochs (|Khalatvan et al1l2008T ). 
Therefore, both SF, magnitude and colour evolution can be 
very complex and is not possible to efficiently probe them. 

We map the evolution of rest-frame galaxy colours up 
to 7-8 billion of years ago, that we determine from our best 
fitted models for 3 values of the look-back time twt corre- 
sponding to redshift z = 0.2, 0.5, 1. In Fig. If 51 we show the 
median values and median deviations for NUV — r and g — r 
of galaxies in the sample. Actually, colour and luminosity 
evolution is stronger for galaxies with a more recent event 
of SF and quenching (i.e. lower t ga i — tAGN)- These galax- 
ies became red only recently, while those with a quenching 
event that occurred many Gyrs ago were already on the red 
sequence. While changes of 2-3 magnitudes in NUV-r are 
observed, g-r changes are <; 1 magnitude, fn addition to our 
reference model, we show the evolution of colours obtained 
using a library of truncated SFs with Z = Zq and r = 3 Gyr 
and two mixed models consisting of a solar single burst and 
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Figure 15. Colour evolution predicted by fit. The different panels refer to values of t ga [ — tAGN accordingly to the same colour coding 
of Figs. ITOl 1111 and 1131 We show the results from our reference model leaving all parameters free to change (star symbols and short- 
dashed error bars), Z = Zq and r = 3Gyr (diamond symbols and continuous error bars), two combined models with a Z = Zq burst 
superimposed to a truncated exponential SF with r = 3 Gyr with a free mass ratio (box symbols and long-dashed error bars) and one 
with single burst amounting to the 75% of the total mass (circle symbols and thin continuous error bars). The black symbols and bars 
are obtained using an unperturbed SF with r and Z free to change (averaged over all galaxies in the sample). At each redshift, the 
different models are artificially shifted to make clearer the appearance. Left panels. NUV-r vs z. Right panels, g-r vs z. 
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a truncated SF with Z — Zq and r = 3 Gyr, leaving free 
the proportion of the two populations, and fixing the sin- 
gle burst to 75% of the total mass. These two latter models 
leave the galaxy colours for t ga i > tAGN almost unchanged, 
while affecting the early history of galaxies producing red- 
der colours, this effect depending on the ratio of the two 
populations. For comparison, we show the median colours 
obtained fitting an unperturbed SF leaving free to change 
tgai, Z and t. 

These predicted tracks show an interesting result con- 
cerning the strength of quenching, that is mainly evident 
in the NUV-r colour. With the exception of galaxies with 
tgai < tAGN and those with an early quenching event, the 
systems with an observed larger colour difference (and thus 
on average large values of t ga i — tAGN) were bluer and more 
star forming at z ~ 0.2. Thus, this result seems to indi- 
cate that in galaxies with a large SF at z ,>0.2, the effect 
of AGN quenching can be stronger, to produce the reddest 
galaxies we observe today. Nonetheless, recent X-ray and op- 
tical selected analysis of high redshift AGNs ([Martin et al.l 
l2007bl . iNandra et aTll200l ISilverman et alJl2008al lb¥ exhibit 
a large fraction of AGN in galaxies at intermediate/blue 
colours, which could be driven by AGN feedback toward 
red-sequence at z ~ 0. 

The magnitudes NUV and r against redshift are plot- 
ted in Fig. 1161 The estimated evolution of the NUV magni- 
tude, being strongly sensitive to RSF, can be different among 
models, and is compared to those obtained using an unper- 
turbed SF. The NUV fluxes are systematically larger (cor- 
responding to lower NUV — r) than those estimated using 
unperturbed SFs. On the contrary, the luminosity in the r 
band does not depend (or at least depends only little) on a 
change in models and the predictions agree quite well. This 
is not surprising, since redder bands are less sensitive to re- 
cent or past SF. 



6 CONCLUSIONS 

In this work, we have attempted to build up a realistic 
model of AGN feedback. We performed a hydrodynami- 
cal simulation to analyze the impact of jets produced by 
AGNs on a inhomogeneous medium, properly representing 
cold gas clouds that form stars in galaxies. In our simula- 
tion, which extends and generalizes a set th at we previously 
performed l|Antonuccio-Delogu fc S ilk 2008), a powerful jet 
Pj = 10 46 ergs" 1 propagates within an inhomogeneous, two- 
phase ISM, containing a realistic distribution of star-forming 
clouds. SFRs in clouds are described using the empirical SK 
prescription, where the SFR depends on the mass density 
of the star forming regions of the clouds. In an early phase, 
the shocks advancing before the expanding cocoon tend to 
compress the cold clouds, without significantly changing the 
temperature of the medium, thus increasing the SF. Later, 
when the cocoon has propagated within the medium, the 
temperature of both the medium and of the clouds increases 
significantly and the mass of clouds is reduced, also due to 
KH instabilities. Thus, at the beginning, a positive feed- 
back increases SF, but the dominant effect is the negative 
feedback that quenches SF in a time of ~ 2 — 3 x 10 6 yr. 
One interesting result of this paper is that for the first time 
a hydrodynamical simulation allows the determination of 



the effect of jets emitted by AGN on SF in galaxies. Pre- 
vious work has insted relied relying on emp irical prescrip- 
tions to take into account AGN feedback (|Granato et al.l 
l200ll . |2004 ICattaneo et~"ai1l2006l . iMartin et alj|2007al ). and 
are thus more robustly supported by our results. 

The high jet power in this simulation is probably the 
main reason of the fact that, at the end of the simulation, 
all the gas within the computational volume is in a physi- 
cal state for which all SF is quenched. For lower injection 
powers, the cocoon expansion will be halted before, and the 
region where SF is quenched will consequently be smaller. 

Based on the results of this simulation, we develop a 
more general model which assumes a SFR for ETGs com- 
posed of a background SF with a more or less extended du- 
ration, that is quenched by a feedback effect such as that an- 
alyzed in the first part of this paper. However, we stress here 
that this general prescription can also describe the effect of 
other quenching mechanisms, like merging, harassment, etc. 
Restricting ourselves to the case of a high power jet, as in the 
present paper, the typical timescales of quenching are very 
small, and no other sources of feedback can stop SF within 
such a short time. Note also that galaxy merging seems to be 
one of the mechanisms which activate AGNs, which after the 
burst induced by t he merger is able to quench the residual 
SF in the remnant ijDi Matteo et alJliooa . iKhalatvan et all 
l2008h . 

Our paradigm is that luminous ETGs are an excep- 
tion within the zoo of galaxies in the universe. We sup- 
pose that the normality in galaxy samples is represented 
by fainter ETGs, that have a protracted SF, but by the ef- 
fect of feedback they are quenched and colours redden. This 
could be one of the missing physical ingredients which could 
explain the short d uration of SF, when exponential SF are 
fitted to data (e.g., lTortora"e t al. 2009). On the contrary, 
semi- analytical galaxy sim ulations predict more protracted 
SFR (|De Lucia et all 120061 ) : supernovae feedback was taken 
into account in these simulations, but the disagreement with 
other observations can be d ue to an improperly a ccounted 
for contribution from AGNs. lSpringel et al] (|2005aT ) simulate 
the effect of AGNs in the merging of two massive and gas- 
rich spiral galaxies, obtaining similar results. They found 
that SF is inhibited with respect to the case without black 
holes; interestingly, without AGNs this kind of merger is not 
able to produce red ETGs, that remain blue (with residual 
SF) even after several Gyrs. On the contrary, the presence 
of black holes redden galaxy colours much faster, giving 
u — r ~ 2.3 in less than 1 Gyr after the beginning of the 
merging process. 

We assume an unperturbed SF law with a scale r = 
1, 3, 5 Gyr, that is quenched at tAGN leaving free to change 
these parameters as well as galaxy ages and metallici- 
ties. We fit these models to observations from a cross- 
matched SDS S +GAL EX catalog. Confirming the results in 
iKavirai et ail (|2007al ). the UV has been shown to be a 
strong indicator of SF, in particular, the UV/optical colour 
NUV — r allows us to select galaxies with different lev- 
els of SF, stopped in the last Gyrs or still in action. We 
show that the quantity t ga i — tAGN is able to describe 
the physical state of galaxies, indicating how much time 
ago SF is stopped. The largest number of galaxies have 
tgai — tAGN ~ 0.5 — 1 Gyr, indicating the necessity of a 
RSF phenomenon (until 1 — 2 Gyr ago), that is quenched 
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Figure 16. Magnitude evolution predicted by fit. See Fig. U5l for details about models and colour or symbol code. Left panels. NUV vs 
z. Right panels, r vs z. 
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by AGNs. Galaxies with larger NUV — r have higher val- 
ues of tgai — tAGN, i.e. SF is quenched early in their history, 
while lower values of NUV — r correspond to galaxies with a 
more recently quenched SF. Finally, galaxies with a SF not 
affected by AGN have the flattest NUV - r colours. This 
distinction is not so clear (or absent) for visible colours. 
These results are shown in Figs. [10] and If f I The epoch of 
the quenching event appears to be correlated with ultra- 
violet flux: galaxies bluer and brighter in NUV band are 
also those with a more recent (or absent) quenching event. 
As shown in Fig. [8] t ga i — tAGN is found to negligibly de- 
pend on the details of background spectral library (e.g., r 
and Z), at least for the bulk of galaxies: thus we obtain a 
robust estimate of the quenching time and RSF. Our anal- 
ysis is less sensitive to the earliest phases of galaxy history, 
notwithstanding the mean evolution of galaxy population 
here analyzed, as shown in Figs. If 5l and ll6l Finally, a short- 
coming of this analysis is related to the timescale of quench- 
ing event, since softer power jets would inhibit SF within a 
larger timescale. Such a slower quenching model would be 
fitted to the colours by a larger t ga i —tAGN, thus our results 
for this parameter have to be interpreted as lower limits. 

One of the most significant results we have found is 
that the NUV-g colour index is very sensitive to the pres- 
ence of a very young stellar population. Typical enhance- 
ment of 2-3 magnitudes of NUV-g are observed, with regard 
to the no-feedback case. Our findings agree with some cur- 
rent simulations that invoke two different modes of AGN 
feedback. A so called 'quasar mode' assumes that during a 
major merger event at high redshift, a fraction of the gas 
accreted by a central black hole i s injected into the gas 
of the host galaxy, qu enching SF (|Springel et al.l l2005al 3, 
|Pi Matteo et al]|2005l ). At later times, another effect is im- 
portant, i.e. the 'radio mode', that is responsible for mak- 
ing galaxies quiescent, and this effect is drive n by low- 
level A GNs (ICroton et al.ll2006l . lBower et al.ll2006l ). More re- 
cently, |Sch^inskrer|ay|2003) speak of a 'truncation mode' 
to indicate AGN feedback at high redshift. At recent epochs, 
no such strong activities or powerful radio jets have been ob- 
served, and so this kind of AGN feedback is often referred 
to as the 'suppression mode'. These aut hors suggest that 
this p rocess could leave a residual SF (e.g. ISchawinski et afl 
2006), but is able to move galaxies along the red sequence. 

Further refinements of our analysis are needed. Firstly, 
it is important to enlarge the sample of galaxies and an- 
alyze a more extended range of magnitudes. The galaxies 
under analysis, due to constraints on magnitude, are rel- 
atively bright with Mb <^ —19. Different kinds of analy- 
sis seem to indicate that brighter and more massive ETGs 
(M B < -20.5 and A/* ^ 10 11 M©) are fundamentally dif- 
ferent from fainter and less massive ones (Mb J> — 20.5 
and M t < 10 11 Mq). In these two different lum inosity and 



mass regimes, the size-luminosity or size-mass ( Shcn ct al 
20031 ) and Faber-Jackson relations (Mat kovic fe Guzman 



20051 ') have differ ent slopes. Also, the Sersi c index is changing 
with luminosity (jPrugniel fc S imicn 199?]), and dark matter 
is shown to have a bivariate behaviour in the two ranges 
(|Tortora et al.l T2009) . suggesting that physical phenomena 
allowing the formation of galaxies of different mass and lu- 
minosity can be variegate. Thus, a future analysis should 
also be directed at studying a wider range of luminosities 
and masses to map the transition from the blue cloud to the 



red sequence. Also to connect our results with the down- 
sizing scenario, we need to enlarge the sample to fainter 
magnitudes, since strong changes in observable quantities 
such as galaxy age, tAGN T and Z are probably visible 
in these luminosity regimes. We are planning to do other 
simulations changing the main input parameters, in order 
to have a more reasonable model of AGN feedback which 
depends, e.g ., on the jet power, lin king it to main galaxy 
observables (|Liu. Jiang, fc Gull2006h . Comparison with sim- 
ulations will allow to quantify how much AGN feedback 
has to be implemented, to make both semi-analytical and 
hydrodynamical simulations able to correctly predict the 
main properties of galaxies we observe. Such a more com- 
plex model of quenching would be parameterized as a func- 
tion of the main parameters of the s ystem and compared 
to other AGN feedback prescriptions ([Granato et alj l200ll . 
|2004 ICattaneo etaHl2006l . IMartin et alj|2007af ). 

Linking observations at low redsh ift, such as those an- 
alyzed here, with h i gh redshift data (IMartin et al.l l2007bl. 
iNandra et ail 12001 [Silverman et al.l l2008al lbh could cer- 
tainly be a powerful way to shed light on the galaxy evo- 
lution scenario, and in particular the role that AGNs have 
in the early and later phases of the SF history of galaxies. 
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